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I. INTRODUCTION 

The potential of static heavy quarks illuminates the most important features of QCD 
dynamics: the asymptotic freedom and confinement. Trying to study subtle electroweak phe- 
nomena in the heavy quark sector of Standard Model, we need quite an accurate quantitative 
understanding of effects caused by the strong interactions. In addition to the perturbative 
calculations for hard contributions, at present there are three general approaches to get a 
systematic description of how the heavy quarks are bound into the hadrons and what are 
the relations between the measured properties of such the hadrons and the characteristics 
of heavy quarks as relevant to the electroweak interactions and QCD. These approaches are 
the Operator Product Expansion (OPE) in the inverse powers of heavy quark mass, the Sum 
Rules (SR) of QCD and the Potential Models (PM) for the systems containing the heavy 
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quarks by exploring various approximations of Bethe-Salpeter equation with the static po- 
tential treated in the framework of effective theory with a power counting in terms of powers 
of the inverse heavy quark mass. The first method is usually exploited in the inclusive esti- 
mates, while the second and third techniques are the frameworks of exclusive calculations. 
The important challenge is a consistency of evaluations obtained in such the ways, that re- 
quires the comparative analysis of calculations. A wide variety of systems and processes for 
the analysis provides a more complete qualitative and quantitative understanding of heavy 
quark dynamics. 

In the leading order of perturbative QCD at short distances and with a linear confining 
term in the infrared region, the potential of static heavy quarks was considered in the Cornell 
model Q, incorporating the simple superposition of both asymptotic limits (the effective 
coulomb and string-like interactions). The observed heavy quarkonia posed in the interme- 
diate distances, where both terms are important for the determination of mass spectra. So, 
the phenomenological approximations of potential (logarithmic one [[| and power law ||), 
taking into account the regularities of such the spectra, were quite successful [§J, while the 
quantities more sensitive to the global properties of potential are the wave functions at the 
origin as related to the leptonic constants and production rates. So, the potentials consistent 
with the asymptotic freedom to one and two loops as well as the linear confinement were 
proposed by Richardson |J, Buchmiiller and Tye ||, respectively. Technically, using a given 
scheme of regularization, say, MS, one has to calculate the perturbative expansion for the 
potential of static quarks. This potential can be written down as the coulomb one with the 
running coupling constant in the so-called V scheme. Thus, the perturbative calculations 
provide us with the matching of MS scheme with V-one. The n loop running of af s re- 
quires the n — 1 loop matching to a v . Note, that initial two coefficients of corresponding 
(3 functions are scheme and gauge independent, while others generally depend. With the 
dynamical fields integrated out, the V scheme is defined in terms of the action depending on 
the static sources (the distance r), so that its (3 function is gauge invariant. The motivation 
by Buchmiiller and Tye was to write down the f3 function of a v consistent with two known 
asymptotic regimes at short and long distances. They proposed the function, which results 
in the effective charge determined by two parameters, only: the perturbative parameter is 
the scale in the running of coupling constant at large virtualities and the nonperturbative 
parameter is the string tension. The necessary inputs are the coefficients of (3 function. 
Two loop results and the one loop matching condition were available to the moment. Re- 
cently, the progress in calculations has provided us with the two loop matching of V and MS 
schemes ||7|||, that can be combined with the three loop running of a™ s . So, the modifica- 
tion of Buchmuller-Tye (BT) potential of static quarks as dictated by the current status of 
perturbative calculations is of great interest. Moreover, to the moment two peculiar ques- 
tions become open. First, the asymptotic perturbative expansion of BT (3 function to the 
third order results in the three loop coefficient, which is wrong even in its sign. Second, the 
elaborated A^g- parameter by BT is in a deep contradiction with the measured value |J. To 
clarify the situation, we are tending to derive the static quark potential consistent with the 
state of the art. 

Thus, our motivation is to combine high order multi-loop calculations of the perturba- 
tive static potential 0,§] with the string tension ansatz. So, we improve the perturbative 
input for the potential model in order to remove the contradiction between the modern high 
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energy data on the QCD coupling constant and the description of heavy quark potential in 
the framework of one-loop Buchmuller-Tye model, which accepts extremely high value of 
coupling constant evolved to the Z mass scale. In other words, if we accept the current nor- 
malization of coupling constant and introduce its value into the Buchmuller-Tye approach 
to the one-loop potential, then we get the contradiction of such the potential with the heavy 
quarkonium mass spectra, certainly, since we find about 200 MeV smaller splitting between 
the 15* and 25 levels in comparison with experimental 580 MeV. This discrepancy cannot be 
removed by the modification of nonperturbative part in the potential with no contradiction 
with the data on the slope of Regge trajectories. Therefore, the modification of perturbative 
input for the model of static potential in QCD is meaningful in this sense even when the 
nonperturbative contribution is conserved in the old string tension form. So, the significant 
improvement of perturbative f3 function for the charge in the coulomb potential is combined 
with the consequent evolution from high virtualities to low ones with taking into account 
the influence of nonperturbative term on the evolution, that becomes essential numerically 
below the scale of 4 GeV. 

We have to emphasize that to the moment of paper by Buchmiiller and Tye || a theory 
for the heavy nonrelativistic QQ pair did not exist. So, the phenomeno logical derivation of 
the static potential including perturbative short- distance and non-perturbative long-distance 
elements made by BT was all one could do. At present, at least for very heavy quarks, such 
a theory does exist in the form of pNRQCD |10| and vNRQCD ||11|| , and we address the 



comparison of static potential model developed in this work with these sound theoretical 
approaches in QCD to the physics of heavy quarkonium. 

Another aspect of this work is devoted to the heavy quark masses. After the potential 
is given, the heavy quark masses incorporated in the corresponding Schrodinger equation 
determine the heavy quarkonium spectra with no ambiguity^. These masses involved in 
the potential model are denoted by tuq. Such the mass should be distinguished from the 
pole mass which is a purely perturbative concept defined unambiguously at each order of 
perturbation theory through the pole of the perturbative heavy quark propagator. Thus, 
we need to test the consistency of estimates for the masses in the QCD potential of static 
quarks and in SR, say. 

In Section 2 we generalize the BT approach to three loops and derive the static potential 
of heavy quarks. The numerical values of potential parameters and their consistency with 
the relevant quantities are considered. The implications for the heavy quark masses, spectra 
of heavy quarkonia and leptonic constants are discussed in Section 3. The obtained results 
are summarized in Conclusion. 



II. QCD AND POTENTIAL OF STATIC QUARKS 

In this section, first, we discuss two regimes for the QCD forces between the static heavy 
quarks: the asymptotic freedom and confinement. Second, we formulate how they can be 



§We deal with the so-called spin-averaged spectra, since the consideration of spin-dependent split- 
ting involves some additional parameters beyond the static potential. 
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combined in a unified (3 function obeyed both limits of small and large QCD couplings. 



A. Perturbative results at short distances 

The static potential is defined in a manifestly gauge invariant way by means of the 



vacuum expectation value of a Wilson loop [12 



V(r) = - lim ^ln(Wr), 
W r = tr7 ? expU0 j> dx^A . (1) 

Here, T is taken as a rectangular loop with time extension T and spatial extension r. The 
gauge fields are path-ordered along the loop, while the color trace is normalized according 
to tr(..) = tr(..)/trl . 

Generally, one introduces the V scheme of QCD coupling constant by the definition of 
QCD potential of static quarks in momentum space as follows: 

Vtf) = -C F ^1, (2) 

while a v can be matched with ajg 

«v(q 2 ) = a m (^) f>„Gu7q 2 ) {^p-J (3) 

oo / ( 2\ \ n 

= ^(q 2 )E«n(^ Li J • (4) 
At present, our knowledge of this expansion^ is restricted by 



a = a = l, a x = — C A - — T F n f , a x = a x + p In — , (5) 

9 9 q 2 



which is the well-known one- loop result, and the recent two-loop calculations which 

gave 

/4343 , „ 7T 4 22 A/ A 9 /1798 56 A/ A „ „, 

02 = (w + ^ - t + t c(3) ) cl - hr + t c(3) ) ^ 

55 le^^^nz+^Tp-n/) , (6) 
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2 2 

S 2 = a 2 + $ In 2 ^ + (Pi + 2p oai ) In (7) 
q 2 q 2 



"On a possible peculiar behaviour in the expansion see ref. [1 
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We have used here the ordinary notations for the SU(iV c ) gauge group: C A = N c , Cp — , v . 
Tp = |. The number of active flavors is denoted by nf. 

After the introduction of a = f-, the (3 function is actually defined by 



din fi 2 



/3(a) = -£/3 n -a"+V), (8) 



n=0 

so that ffl tl = ffiff and 

^ = pf- ai pf s + (a 2 -al)PT s (9) 
,618 + 24203) ll(167T 2 -7r 4 ) ^ 3 
~ V 9 12 / 

_ ( 445 + 704 C (3) + lg^ )clrfn/ 

+i± 224C(3) C4(rFn/)2 _ ^-^CP) ^^ 

+2Cf7> n . / + 184 ~ V^- C^Tpn,) 2 . (10) 



The coefficients of (3 function, calculated in the MS scheme |13] are given by 

PT= l ^C A - A -T F n f , (11) 

34 Of) 

0ms = _ C 2 a _ ACFTFUf _ -C A T F n f , (12) 

ft* = + 2C F Tpn f - 2 -fc A C F Tpn f - ^C\T F n s 

+^Cp(T F n f y + l ^-C A {T F n f f. (13) 
The Fourier transform results in the position-space potential J7J 

Qm s(^ 2 ) L «ms(^ 2 
r \ Air 

,2\ n , ^2 



y( r ) = i + ^ii(2/3 ln(/irO + <■ 

r \ An v 

+ (^T 1 ) 2 (^(41n 2 (^) + y) (14) 

+2(/3 a + 2#)Oi) ln(//r') + a 2 ) + . . . J 

with r' = r exp^p). Defining the new running coupling constant, depending on the distance, 

V{r) = -Cp^^l. (15) 
r 

we can calculate its /^-function from fll4|), so that |7j] 

& v = & v + y/3 3 , (16) 

and the minor coefficients /3q 1 are equal to the scheme- independent values given above. 
To normalize the couplings, we use (§) at q 2 = m 2 z . 
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B. Confining term 

The nonperturbative behaviour of QCD forces between the static heavy quarks at long 
distances r is usually represented by the linear potential (see discussion in ref. | fH| ) 

V con \r) =k-r, (17) 

which corresponds to the square-law limit for the Wilson loop. 

We can represent this potential in terms of constant chromoelectric field between the 
sources posed in the fundamental representation of SU(iV c ). So, in the Fock-Schwinger 
gauge of fixed point 

x M • A"{x) = 0, 



we can represent the gluon field by means of strength tensor [IS 



A^x) « -l x u G utl (0), 

so that for the static quarks separated by the distance r 

0,(0) G- (0) Q 3 {0) = r -yET? p 

where the heavy quark fields are normalized to unit. Then, the confining potential is written 
down as 

V c ™\r) = l -g s C F E-r. 

Supposing, that the same strength of the field is responsible for the forming of gluon con- 
densate, and introducing the colored sources rii, which have to be averaged in the vacuum, 
we can easily find [[Uj 

(GW = -4(G^ (0)G^ (0)) =AC F E 2 (nn) J 

where we have supposed the relation 

(nT a T b n) = - (nT a n ■ nT b n) , (18) 

which ensures that the sources conserve the massless of the gluon, and, hence, the gauge 
invariancep]. Further, it is evident that 

(nT a T b n)=C F -f—(nn). 

c 

Then, we conclude that the relation between the strength E and the string tension depends 
on the normalization of vacuum sources n«. We put 



(mrij) = m5. 



t^The mass term generated by the sources should be equal to C ~ A^A b v [nT a T b n + nT a n ■ nT b n], 
so that the averaging of sources yields zero, if we suppose (|i"8|). 
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where n\ denotes the number of light stochastic flavors, which is the free parameter of such 
the representation. Of course, the value of n\ should be finite even in the case of pure gluo- 
dynamics with no light quarks in the infrared region. Moreover, the light quark loops could 
cause the breaking of gluon string, i.e. the strong decays of higher excitations. We assume 
that ni is basically determined by the gluon dynamics (i.e. the number of colors), and it 
slightly correlates with the number of quark flavors. After a simple consideration of poten- 
tial strength between two colored sources in the fundamental and adjoint representations, 
i.e. the color factors in front of single gluon coulomb potential, we assume that in the pure 
gluodynamics the number of stochastic sources substituting for the vacuum gluons can be 
accepted in the formP] 

1 C A 3 1 
ni= N- c C- F = l = l nU 
where the factor 1/N C normalizes the source to unit, and Ca/Cf is the appropriate ratio of 
color charges. To the moment, the shift of n; in QCD with light quarks is not explicitly fixed, 
while the lattice calculations shown that the dependence of string tension on the number of 
light quarks is weak |17| . Finally, we find for the linear term of the potential 



71 

VCF~N c ni 



I = — C Fy j{^GU) = ^= C FS j(^Gl v ). (19) 



The k term is usually represented through a parameter a' P as 

k 1 



27T«p 



Buchmiiller and Tye put a' P = 1.04 GeV 2 , which we use throughout of this paper. This 
value of tension, that is related with a slope of Regge trajectories, can be compared with 
the estimate following from (0). At (^G%,) = (1.6 ± 0.1) ■ 10~ 2 GeV 4 @ we have found 

a' P = 1.04 ±0.03 GeV" 2 , 

which is in a good agreement with the fixed valueQ. 

The form of flT7|) corresponds to the limit, when at low virtualities q 2 — > the coupling 
QV tends to 



«v(q 2 ) -> 



K 
q 2 



so that 



lnq 2 



-«v(q 2 ), (20) 



"This assumption corresponds to the definition of vacuum properties in QCD in terms of notations 
under the consideration, which is in agreement with the value of gluon condensate and Regge 
trajectories slope. 

^The ambiguity in the choice of n\ can change the appropriate value of gluon condensate. 
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which gives the confinement asymptotics for the (3 V function. 

A special comment should be done on the role of linear term in the potential. Considering 
the power corrections, which can be attributed to various sources such as the renormalon, 
topological effects caused by monopoles and vorteces, deviations from the operator product 
expansion, the authors of ]19j argued that this term responsible for the quark confinement 
can contribute at short distances, too. This conclusion is essentially different from the point 
of view based on the notion about a low energy phase transition leading to the condensation 
of gluons and quarks. This condensation provides the forming of chromoelectric string 
between the static quarks. Thus, at short distances (or high virtualities q 2 ) one could expect 
the decomposition of condensates, that means the scale of confinement (or the string tension) 
should disappear from the physical quantities at large q 2 . In contrast, the nonperturbative 
scale can contribute as the factor in front of power corrections 1/q 2 even at q 2 — > oo. So, in 
[1^1 several indications of linear term contribution at small distances were considered. We 



repeat the items relevant to the question on the static potential here. 



First, the lattice simulation [fJUJ does not show any change in the slope of the full QQ 
potential as the distances are changed from the largest to the smallest ones where the 
coulombic part becomes dominant. Hence, no rapid energy jump, characteristic for the phase 
transition, is found on the lattice. An explicit subtraction of the perturbative corrections at 



small distances from the potential in the lattice gluodynamics was performed in [21]. This 
procedure gives an essential nonzero linear term at very small distances. 



Second, there are the lattice measurements |22] of the fine splitting in the heavy quarko 



nium levels as a function of the heavy quark mass. The approach by Voloshin and Leutwyler 
|16|1 predicts a particular pattern of such the dependence. Indeed, the multipole expansion 



of heavy quarkonium interaction with the external gluon field leads to the dominant contri- 
bution by the second order of chromoelectric dipole. Therefore, the quark distance squared 
appears as the leading term in the perturbation due to soft gluons at short distances. These 
predictions are very different from the evaluations based on the static quark potential with 
the linear term, such as the potential by Buchmiiller and Tye 0. The numerical results 
from the lattice simulations favor the linear correction to the potential at short distances. 

Third, an interesting manifestation of short strings might be the power corrections to 
current correlation functions IL(g 2 ). Calculations of a relevant coefficient in front of the 1/q 2 
terms involve the model assumptions. So, in |23[ it was suggested to simulate this power 



correction by a tachyonic gluon mass. The tachyonic mass can imitate the stringy piece 
in the potential at short distances. Rather unexpectedly, the use of the tachyonic gluon 
mass (m 2 = —0.5 GeV 2 ) explains well the behavior of Hj(q 2 ) in various channels. This fact 
implies again we see the confirmation of short distance linear term in the potential. 

Thus, we do not involve any additional assumptions on the possible scale and properties 
of quark-gluon condensate decomposition at short distances in the description of static 
potential in QCD. 



C. Unified j3 function and potential 

Buchmiiller and Tye supposed the following procedure for the reconstruction of f3 function 
in the whole region of charge variation by the known limits of asymptotic freedom to a given 



S 



order in a s and confinement regime. So, in the framework of asymptotic perturbative theory 
(PT) to one loop, the At is transformed to the Richardson one, 



At (a) 



A<* 2 Aiich(a) ^ oa 2 (i 

— exp 



1 

'/3oo 



(21) 



The Richardson function has the essential peculiarity at o — > 0, so that the expansion is the 
asymptotic series in a. At a — ► oo the f3 function tends to the confinement limit represented 
in©- 

To the two loop accuracy, following in the same way results in the j3 function by 
Buchmiiller-Tye, 



1 + A 



1 



1 



Ar(a) p a 2 ffia /3 BT (a) /5 a 2 fl-exp 



i 

"ft) a 



+ §z^M~H- (22) 



The exponential factor in the second term contributes to the next order in a at small a, so 
that the perturbative limit is restored. However, we can easily find that third coefficient of 
/?bt function is equal to 

/3 2> BT = ^(/?l-^0), 

Po 

and it is negative at the chosen value of I = 24 || , which is in contradiction with the recent 
result [0|§, shown in ([10]). 

To incorporate the three loop results on the perturbative (3 function, we introduce 



1 



At (a) 
1 



[3 a 



1 A + 
+ 



ft) 



A 2 a 



/3(a) (3 Q a 2 (l - exp 



i 

'ft a 



+ 



A + 



ft) 



A 2 a 



exp 



l 2 a 2 



(23) 



where again the exponential factor in the second term contributes to the next order in a —>■ 0. 
In the perturbative limit the usual solution 



a(/i 2 



A) lug 



A i , , ^ 2 , A 2 

A 2 ln^ 



A 2 ' A 4 In 2 £ 



< ^2 ^2 

In In — — In In — 



A 2 



A 2 



1 + 



A V A ^ 
A 2 , 



(24) 



is valid. Using the asymptotic limits of and (0), one can get the equations for any (3 
function, satisfying these boundary conditions, as follows: 



In 



l A f a ^> 

— — ^- + -2 lnAa(/^ 2 ) + / dx 
Aa(/i 2 ) A -A) 

In — = In o(/i 2 ) + / 



A 2 
A" 

7 2 



Pi + i 



A^ 2 Po x Pi x ) 



dx 



1 1 

a; /3(a;) 



(25) 
(26) 



In general, at a given (3 function, Eqs.(^) and ( [26]) determine the connection between the 
scale A and the parameter of linear potential K, 



k = 2nC F K. 
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Supposing (|23|) we can easily integrate out ( j25|) to get the implicit solution of charge depen- 
dence on the scale 



In 4 = In 



A 2 



cxp 



Joa(^ 



- 1 



+ 



In 



/3o\/2 1 



1e + Ei 



Z 2 aV) 



+ 



1 V_2 



Erf 



v/2 



(27) 



where Ei[x] = cit t — 1 exp[— t] is the exponential integral, and Erf[x] = ^= exp[— t 2 ] 
is the error function. 

Eq.(|27|) can be inverted by the iteration procedure as it was explored in the derivation 
of (pip. So, the approximate solution of ( |2"7D has the following form: 

a(/i 2 ) ' 



A) In (l+r^ 2 )^) 



(28) 



where 



I w 



cxp 



IE +Ei 



/3J/3 — /? 2 \/ 2 



Erf 



I ai(/i 2 



while di is obtained in two iterations 



fa In (i + Vi (jjl% 



*7l(/* 2 ) 
Oo(^ 2 y 



exp 



/3 



2/? 2 



^ 7s + Ei 



/3 ln(l + ^)' 



(29) 

(30) 

(31) 
(32) 



Taking the limit of /i 2 — > we find the relation 



hxA<K z C F a' P k* = ln/5 + 7s + 



2/3 2 



/TO -Pi \2 

(31 I 



(33) 



which completely fixes the f3 function and charge in terms of scale A and the slope a' P , since 
we have expressed the parameter Z in terms of above quantities. 

Remember, that at /i 2 — > 00 the perturbative expression (|54]) becomes valid as the limit 
of effective charge (p8|). 

To the moment we are ready to discuss the numerical values of parameters. 



D. Setting the scales 

As we have already mentioned the slope of Regge trajectories, determining the linear 
part of potential, is fixed as 

a' P = 1.04 GeV~ 2 . 
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We use also the measured value of QCD coupling constant || and pose 

af (m|) = 0.123, 

as the basic input of the potential. 

At the given choice of normalization value for the QCD coupling constant we get the 
scale A^_ 5 ~ 273 MeV, which certainly differs from the world average value resulted in 
the analysis of PDG @, where A^ =5 ~ 208^23 MeV, that corresponds to the coupling 
constant af s (m 2 z ) = 0.1181 ± 0.002 Q. However, this average value including various data 
is generally determined by the most precise measurements: the data on the hadronic events 
in the peak of Z boson at LEP (the hadronic width), the decays of r lepton, the data on 
the deep inelastic scattering (DIS) for leptons off nucleons and the lattice simulations for 
the systems of heavy quarkonia. In this set of estimates, the high energy measurements at 
LEP for Z and at HERA for the evolution of nucleon structure functions give the average 
values ctf^(mi) = 0.123 ± 0.004 and cr^(m|) = 0.122 ± 0.004, respectively, while the 
evolution of structure functions at low virtualities, where an ambiguity in the description 
of nonperturbative effects and contributions of higher twists are essential, as well as the 
energy-dependent sum rules for the structure functions at low energies significantly displace 
down the common average value for the coupling constant extracted from the DIS data. 
Thus, we argue that the methodical uncertainty for such the averaging is underestimated, 
since the low-energy data have got some little calculated sources of theoretical uncertainties. 
The analysis of data on the decays of r lepton resulting in a™ s (m 2 z ) = 0.121 ± 0.003, is 
based on the sum rules, where the control of nonperturbative corrections is much better 
than in DIS, though there are some theoretical problems on the formulation of sum rules 
in the region of physical states in contrast with the classic variant of sum rules in the deep 
euclidean region. Finally, the lattice simulations investigate the splitting between the states 
of heavy quarkonia, i.e. they operate with the low-energy data and rely on the approximation 
with the zero number of light quarks = or nj = 2 under the extrapolation to both 
the real number of nj = 3 and the region of high virtualities due to the evolution. A 
high accuracy of such lattice estimates is announced. As we have seen the spectroscopic 
characteristics for the systems of heavy quarks need an extremely careful interpretation, since 
the evolution of potential parameters from the region of bound states to the high virtualities 
is affected by the nonperturbative factors. Thus, we see that the normalization value of 
QCD coupling constant accepted above agrees with the direct high-energy measurements, 
while the data obtained at low energies allow the agreement, if we take into account their 
systematic uncertainties, which are not well estimated. 

Note that the decrease of normalization value to af s (m z ) = 0.120, for example, leads 
to the discrepancy with the data on the splitting of heavy quarkonium masses between the 
levels of 15* and 25* states, which is very sensitive to the normalization of QCD coupling 
constant, so that instead of M{2S) — M(1S) ~ 580 MeV we get the value which is less by 
about 100 MeV. In this respect, the variation of other dimensional parameter, the Regge 
trajectory slope, from the accepted value of a' P = 1.04 GeV" 2 to a' P = 0.87 GeV -2 leads 
to unessential change in both the splitting and the corresponding value for the scale in the 
coupling constant evolved to low virtualities. 

Then, we evaluate 

a v (m|) « 0.1306, 
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and put it as the normalization point for a(m|) = a v (m^)/(47r). Further, we find the 
following values of A for the effective charge, depending on the number of active flavors: 



A n/=3 = 643.48 MeV, I = 56, (34) 
A n/=4 = 495.24 MeV, I = 37.876, (35) 
A n/=5 = 369.99 MeV, I = 23.8967, (36) 

where we set the threshold values for the switching the number of flavors to be equal to 
m 5 = 4.6 GeV and = 1.5 GeV. After such the fixing the momentum space dependence 
of the charge, we perform the Fourier transform to get 

8C 

V(r) = k ■ r u(r), (37) 



M ( r ) = /_ — ( a (<? 2 ) - -j J sin (9 



with 

poo 

u(r) = I , . 

q \ q\ 

which is calculated numerically at r > 0.01 fm and represented in the MATHEMATICA file 
in the format of notebook at the site |http:/ /www. ihep.su/~kiselev/Potentia I. nb| . 

Note, that at short distances the potential behaviour is purely perturbative, so that at 
r < 0.01 fm we put 

nr) = _ Cp M2M (38) 

where the running a v (l/r 2 ) is given by eq.(^) with the appropriate value of /3J at rif = 5, 
and with the matching with the potential (|37|) at r s = 0.01 fm, where we have found 



a v (l/r s 2 ) = 0.22213, 

which implies K f=5 = 617.42 MeV. 

Thus, we have completely determined the model for the potential of static heavy quarks 
in QCD. In Fig. [I] we present it versus the distance between the quarks. As we can see the 
potential is very close to what was obtained in the Cornell model in the phenomenological 
manner by fitting the mass spectra of heavy quarkonia. 

The visual deviation between the QCD potential derived and the Cornell model at long 
distances is caused by a numerical difference in the choice of string tension: we adopt the 
value given by Buchmiiller and Tye, while in the Cornell model the tension is slightly greater 
than that of we have used. A more essential point is the deviation between the potentials 
at short distances (see Fig. ||), because of clear physical reason, the running of coupling 
constant in QCD in contrast to the constant effective value in the Cornell model. 
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FIG. 1. The potential of static heavy quarks in QCD (solid line) in comparison with the Cornell 
model (dashed line) (up to an additive shift of energy scale). 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 



r, fm 



FIG. 2. The potential of static heavy quarks in QCD (solid line) in comparison with the Cornell 
model (dashed line, up to an additive shift of energy scale) and the difference between them (upper 
curve) at short distances as caused by the running of coupling in QCD. 
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FIG. 3. The differences between the (3 functions versus the effective charge. The value of 
{(3 — /?bt)/|/3| is shown by the solid line, ((3 — /3Ri c h)/|/3| is given by the short-dashed line, and 
(/3bt — AiichVI/^BTl is represented by the long-dashed line. 



To compare, we show the differences between the f3 functions (21), ( p2"D and (23) in Fig. 



|3] at the fixed values of I and n/ = 3. Wee see that the asymptotic perturbative expansion of 
[3 at a — > dominates at a < do, where do ~ 0.03 corresponding to a v ,o ~ 0.37. This value 
of coupling a Vi o coincides with the effective coulomb constant used in the Cornell model. At 
larger values of coupling the contributions related with the confinement regime are essential. 

Two comments are to the point. First, the resulting potential is obtained by the per- 
turbative normalization to the measured value of af s (m 2 z ) as combined with the three- loop 
evolution to the lower virtualities. Second, the running of coupling constant is modified 
(numerically the deviation from the perturbative regime begins at \x < 3 — 4 GeV) to reach 
the confinement limit at /i — > 0, so that the perturbative connection between the scales A 
and A MS is broken at virtualities under touch by the charmed and bottom quarks, that was 
the reason for the error in the assignment of A MS by Buchmuller and Tye. 



III. HEAVY QUARK MASSES AND LEPTONIC CONSTANTS 

Considering the characteristics of heavy quark bound states we should emphasize a sig- 
nificant necessity to certainly separate two distinct theoretical problems. The first problem 
is the calculation of heavy quark potential, where the leading approximation is the static 
limit of rriQ — > oo in the operator product expansion over the powers of inverse heavy quark 
mass, we have considered this problem in the previous section. The other problem is the 
calculation of bound state masses. In the heavy quarkonium the kinetic energy of quark 
motion is comparable with the potential energy. So, the leading approximation for the ef- 
fective lagrangian in the operator product expansion over the inverse heavy quark mass is 
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the sum of nonrelativistic kinetic term and the static potential, which give the dominant 
contribution in the Schrodinger equation for the bound states. Corrections are relativistic 
terms in the kinetic energy and perturbations of the static potential in the form of operators 
suppressed by the inverse powers of heavy quark mass, as well as nonpotential retardation 
effects. The magnitude of such the corrections can be restricted numerically, that leads to 
a systematic uncertainty in the calculations of mass spectra for the heavy quarkonia in the 
framework of potential approach with the static potential. 



A. Masses 

The determination of potential provides us with the extraction of heavy quark masses 
in the static approximation by comparison of heavy quarkonium mass-spectra with the 
calculated ones. The predicted charmonium and bottomonium masses are presented in 
Tablesf**| | and |T1] at the following values of heavy quark masses in the potential approach 

m v c = 1.468 GeV, m Y b = 4.873 GeV, (39) 

with no taking into account relativistic corrections, which can be sizable for the charmonium 
(say, AM(cc) ~ 40 MeV). To the moment, the only measured splitting of nS'-levels is that 
of rj c and J/ip, which allows us to evaluate the so-called spin- averaged mass 

M(1S) = (3Mj^ + M Vc )/A. 

Supposing the simple relation [[|: M{ns) = Mv(nS) — j-{Mj/^ — M Vc ), we estimate also 
the expected values for the excited states with an accuracy better than 10 MeV, we think. 
For the P-wave levels we explore the masses 

M(P) = M x + i(M 2 - M ) + ^(M 2 -M x + 2(M - Mi)), 

where Mj denotes the mass of state with the total spin J and the sum of quark spins S — 1, 
and we have supposed the spin-dependent forces in the form 

V SD = A(L ■ S) + B(L ■ S) 2 - ii?L 2 ■ S 2 , 

where the third term corresponds to the third term in the above expression for M(P) 



TABLE I. The masses of charmonium as predicted in the present paper (K 2 0) in comparison 
with the experimental data elaborated as described in the text. 



state, nL 


M, K 2 


M, exp. 


state, nL 


M, K 2 


M, exp. 


IS 


3.068 


3.068 


2P 


3.493 


3.525 


2S 


3.670 


3.671 


3P 


3.941 




3S 


4.092 


4.040 


3D 


3.785 


3.770 



***We suppose that the V'(3770)-state is a mixture of 3S and 3D levels with unessential shift of 
3Z?-mass. 
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and it results in the L-dependent shift of levels. 
We have supposed also 



a s {m b ) m 2 c \R b ^ 
a s (m c ) ml |i2gc(0)| : 



M T - M Vb w np^^|^riL( Mj/V) - M Vc ) « 56 MeV. 



We have found that the sizes of quarkonia are the same as they were predicted by 
Buchmiiller and Tye, while the masses of states are slightly different since we have used the 
other prescription for the input values of ground state masses: 

Msc(lS) = 3.068 GeV, M bb (lS) = 9.446 GeV. 

Then, we predict the masses of be q uarkoniumf^l, as shown in Table [TTI|. The calcu- 
lated values of masses agree with those of estimated in the Buchmiiller-Tye and Martin 
potentials ||25|| . The wave functions at the origin are related with the production rates of 
heavy quarkonia. These parameters are close to what was predicted in the BT potential, but 
slightly smaller because of both the change in the charmed quark mass and the asymptotic 
behaviour at r — > 0. 

TABLE II. The masses of bottomonium as predicted in the present paper (K 2 0) in comparison 
with the experimental data elaborated as described in the text. 



state, nL 


M, K 2 


M, exp. 


state, nL 


M, K 2 


M, exp. 


IS 


9.446 


9.446 


2P 


9.879 


9.900 


2S 


10.004 


10.013 


3P 


10.239 


10.260 


3S 


10.340 


10.348 


3D 


10.132 




4S 


10.606 


10.575 


5S 


10.835 


10.865 



TABLE III. The masses of be as predicted in the present paper (K 2 0) in comparison with the 
experimental data. 



state, nL 


M, K 2 


M, exp. 


state, nL 


M, K 2 


M, exp. 


IS 


6.322 


6.40 


2P 


6.739 




2S 


6.895 




3P 


7.148 




3S 


7.279 




3D 


7.013 





tt^The experimental error in the ground state mass is still large, 6M = ±0.39 GeV 
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To the moment we have fixed the potential masses of heavy quarks (|39|) as independent 
of scale. To compare with the masses evaluated in the framework of QCD sum rules, we note 
that in the sum rules for the heavy quarkonia one usually explores the NRQCD [f26] with 
the perturbative potential (|14]) explicitly dependent of the normalization point /i (referred 
as yLi S oft in the SR). We have checked that at short distances and high /i so f t the perturbative 
potential (|14]) and that of present paper coincide with each other, while a deviation appears 
at r > l//i so f t . However, at the distances characteristic for the ground states of heavy 
quarkonia: (r^ns)) ~ 0.22 fm and (r 5c (is)) ~ 0.42 fm, the shape of the potential can be 
approximated by the perturbative term at /i soft = 1.5 — 2.0 GeV (see Figs. [| and ^|) with 
the additive shift of energy scale SV(fi so f t ), which is defined by the expression 

5V(fi soit ) = [V(r) - V peH (r; /i soft )] r= , (40) 

where the parameter ( has been put in the region of £ = 1 — 2, where the energy shift 
5V has got a little variation about 30-40 MeV, which is, on the first hand, a characteristic 
uncertainty of potential approach, and on the other hand, it points to a similar form of 
perturbative potential with the calculated model potential in the region of distance variation. 
The dependence of energy shift is represented in Fig. [| 
So, if we redefine the heavy quark masses^ by 

m pole H, c = < c + Isvfji), 

the solution of Schrodinger equation with the perturbative potential and m pole (/z) results in 
the quarkonia masses close to the experimental values. Thus, we have matched the values of 
potential masses m v in the QCD potential with the perturbative pole masses standing in the 
two-loop calculations. We stress that the dependence on the soft scale in both the energy 
shift 5V(fi) and the pole mass m pole (/i) does not reflect a nonzero anomalous dimension, since 
these quantities are renormalization group invariants. This scale dependence is due to the 
truncation of perturbative expansion, wherein the coefficients in front of powers of coupling 
constant can contain the factorial growth (the renormalon), so that even at the scale close to 
the charmed quark mass the infrared singularity in the running coupling constant of QCD 
provides the significant custodial scale dependence. 

Numerically, we estimate the running masses m{m) in the MS scheme using the two and 
three-loop relations with the pole mass derived in |29| , |30|| and adjusting the scale /z so f t to be 
equal to rh. So, in two loops P§) we get 



m c (m c ) 2loovs = 1.40 ± 0.09 GeV, m b (m fe ) 2l00ps = 4.20 ± 0.06 GeV, 



■^This redefinition is the indication of perturbative renormalon (see review in [27]). Indeed, there 
are two sources for the deviation 5V. The first is the linear confining term in the potential of static 
quarks. However, it is a small fraction of 5V. The second source is the infrared singularity in the 
perturbative running coupling. One can easily find that subtracting the singular term of the form 
~ — ^-r from 5V results in a small value slowly depending on u so f t . In the effective theory for the 
nonrelativistic heavy quarks, the subtraction that connects the pole mass and the threshold mass 
can be calculated explicitly (see [^] and references therein). 
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FIG. 4. The potential of static heavy quarks in QCD (solid line) in comparison with the pertur- 
bative term (|i"4|) at fi = 1.5 GeV (dashed line) (up to an additive shift of energy scale). 
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FIG. 5. The potential of static heavy quarks in QCD (solid line) in comparison with the pertur- 
bative term (|l~4|) at fi = 2.0 GeV (dashed line) (up to an additive shift of energy scale). 
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FIG. 6. The value of additive shift of energy scale to match the perturbative /U-dependent po- 
tential with that of calculated in QCD. The solid and dashed lines correspond to the two and one 
loop matching. The points give the result of sum rules for bottomonium in comparison with the 
dotted curve following from the relation between the running and pole masses at scale fi. 



while the three-loop approximation |K| , which is consistent with the three- loop evolution of 
coupling constant, results in slightly smaller masses, especially, for the charmed quark, where 
the uncertainty of estimate increases because of stronger sensitivity of quantities involved 
to the scale variation, 



m r [m 



C J 3 loops 



1.17 ± 0.10 GeV, rh b {m b ) 3l00ps = 4.15 ± 0.06 GeV, 



which are in agreement with the various estimates in the sum rules on m& [|3~l~H34jl and m. 
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The uncertainty of estimates is determined by the deviations in the calculations of heavy 
quarkonium masses \8M = 20 MeV (as shown in Tables | and ||) and the error in the 
extraction of SV mentioned above. The uncertainty in the running mass of charmed quark 
is slightly larger than in the bottom mass, since, in addition, its value is more sensitive to a 
small variation of scale, pole mass and energy shift. 

Note, that the calculations in the framework of sum rules were performed for the 6-quark 
mass in both the full QCD |36|] and the effective theory of nonrelativistic heavy quarks 
NRQCD |pi- |5^J57| . The mass extraction of Ref. |57|] has been carried out in the nonrel- 
ativistic effective theory at next-to-leading order (NLO), whereas Refs. |31|-p4| carried out 



^Note, there is the difference between the usually quoted values of m(m) and m(m 



pole\ 
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NNLO analyses in the same framework. The calculations in the nonrelativistic effective 
theory are the calculations in the framework of first principles in QCD, where the results of 



full QCD are determined in a systematic expansion in a s and the velocity. In Ref. |36| the 
analysis has also been carried out in the nonrelativistic situation, but no systematic expan- 
sion in a s and the velocity has been carried out. That the results for the MS mass obtained 



in [36| agree with the other analyses is not understood and requires further examination (see 



so 



the conclusions of j3"B"||). 

Recently, the charmed quark mass was evaluated from the NRQCD sum rules in [|39 
that the result on the running mass is in a good agreement with the value given above, too. 
There is a recent sum rule extraction [[!(]] of the bottom MS mass, where the charmed quark 



mass effects are also included. The estimate of potential approach under consideration is in 
a good agreement with this recent SR result. 



In |}3f the dependence of 'pole' mass on the scale /i so f t was explicitly calculated in the 
N 2 LO. The uncertainty of mass extraction from the sum rules for bottomonium was given by 
0.1 GeV for the running MS mass and 0.06 GeV for the low-energy running mass ('kinetic' 
mass). The result on the 6-quark pole mass depends on both the scale of calculations and 
the order in a s of perturbative QCD. To compare the results in the sum rules with those of 
given in the present paper we fix the order in a s by the two-loop corrections. Then we have 
found that, say, at /A so f t = 2.5 GeV the results of estimates in the perturbative potential 
approach and in the framework of sum rules are the same within the uncertainty mentioned. 
So, putting the above value as the matching point we show the sum rule results in the form 
of energy shift in Fig. |6|. For the sake of representability in Fig. ^| we show the //.-dependent 
'pole' mass extracted in [[33] with the uncertainty of 5m = 80 MeV, which is characteristic 
inherent error for the short- distance masses in the analysis of [[33|]. Despite of various choice 
for the normalization of QCD coupling constant (in af s (m 2 z ) = 0.118), we see a good 
agreement between the //-dependencies of both the energy shift in the perturbative potential 
with respect to the static potential of QCD and the variation of perturbative 'pole' mass of 
6-quark in the sum rules of QCD. As for the one-loop matching of perturbative potential, 
we mention only that the corresponding sum rules in the NLO give the value of energy shift 
close to zero at /i so f t > 2 GeV within the uncertainty of the method, and this estimate is 
consistent with the result of potential approach as shown in Fig. §. Thus, the energy shift of 
perturbative potential with the two-loop matching of V and MS schemes indicates the form 
of QCD potential in agreement with the corresponding soft scale dependence of perturbative 
pole mass in sum rules of QCD for the bottomonium. 

To the moment we can compare the obtained /^-dependence of 'pole' mass with the 



relation between the running MS-mass of heavy quark and the pole mass derived in |ET 
where we find 

m po,e = m(/i) | 1 + Cl ( /i )^!^ + cM I ^L-^-L ] I ill) 



47T 




with 



Cl (/i)=C F (4 + 3L), (42) 
c 2 (/i) = C F C A - 8C(2) - 4/ 8 (l) + ^f-L + y/ 2 
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-C F T F n f + 8C(2) + y L + 2L 2 ) (43) 

27 Q \ 

+CJ. + 30C(2) + 8/ 3 (l) + -L + -L 2 j - 12C F I>(1 - 2((2)), 

where J 3 (l) = |C(3) - 6C(2)ln2, and L = 2\n(/j/m pole ). At /i = m pole , the result of 
29fl is reproduced. We check that the logs in the definitions of can be removed by 
the expression of running values m(u) and a^ s (/i) in terms of m(m pole ) and af s (m pole ) 
in (|4T|) . Nevertheless, we find that the explicit /i-dependence in (|4i~|) repeats the form of 
renormalon contribution as we see it in the perturbative potential, where the similar effect 
takes place because of both the truncation of perturbative series and the infrared pole in 
the running coupling constant of QCD. Following (£|l]), we show the value of difference 
2(m pole (/i) — nib) in Fig. || at m{m) = 4.3 GeV. We see that, first, the results of QCD sum 
rules in [|3|] agree with the values expected from (|41~1) , and second, the /i-dependent shift 
of pole mass approximately coincides with the shift of perturbative potential with respect 
to the static QCD potential free off renormalon ambiguity caused by infrared singularity of 
perturbative coupling constant at finite energy scale. This fact implies the cancellation of 
infrared uncertainties. Thus, we can define the unambiguous pole mass by 



^pole = m pole 



1 



(Ai) - -6V(im), (44) 



2 



where we use the pole mass of (fy]). The basis for the validity of ( f44|) was observed in [42 



where in the context of perturbative bottom mass extractions, the cancellation of the leading 
renormalon at u = 1/2 of the Borel plane in the total static perturbative energy of a heavy 
QQ pair was shown. 

We find that for the bottom quark the defined mass is given by the value of mass extracted 
from the potential approach 



~ pole _ v 

ml ps m h , 



with the accuracy about 80 MeV. 



B. Heavy quark masses and pNRQCD 

In this section we discuss the modern development in the theory of heavy quarkonium 
QQ' on the basis of effective theory called pNRQCD [|TUJ, naturally incorporating the po- 



tential interactions between the heavy quarks and external ultrasoft fields in QCD, and 
compare the pNRQCD results with the values of heavy quark masses obtained above in the 
QCD potential of static quarks. 

First, pNRQCD argues that in the heavy quarkonium the nonrelativistic motion of heavy 
quarks inside the bound state allows us to introduce three actual physical scales: the heavy 
quark mass m, the soft scale of heavy quark momentum inside the hadron mv and ultrasoft 
scale of energy mv 2 , which are distinctly separated by a small parameter v being the velocity 
of heavy quark. After the matching with full QCD at a hard scale /^hard ~ in, in NRQCD 
the hard fields are integrated out, that results in the perturbative Wilson coefficients of OPE 
in the effective theory, and we deal with the heavy quarks interacting with the gluons at 
virtualities /!f ac t,soft about mv. In order to consider the heavy quark fields at lower fi up 
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to mv 2 we should introduce the effective lagrangian of pNRQCD, where the soft fields are 
integrated out, and we deal with the potential interaction of heavy quarks and the ultrasoft 
external gluon fields in the framework of multipole expansion. The matching of pNRQCD 
with NRQCD takes place at a scale /if act ~ mv. Recently, the effective theory of vNRQCD 



was formulated in |TTJ, using the velocity renormalization group |g3| to match the vNRQCD 
operators with the full QCD at a scale about m with the single-step evolution to a soft 
scale, which can be either mv or mv 2 . The current status of vNRQCD provides us with 
the one-loop matching of heavy quark potential to order v 2 , i.e. up to spin- dependent 1/m 2 
terms, which are beyond the current consideration. Therefore, we concentrate our discussion 
on pNRQCD. 

The pNRQCD lagrangian has the following form: 

p2 p 2 4 ^(1) y{2) \ 

£ p nrqcd = Tr ■! S T I id — — + -=— , — V s {r) - s — + . . . S 

I \ \ i > i i»i ! Mi'' <yy~y Tfi i 




-V (r --2 °— + ... O (45) 

mm 1 J J 

V B {r), 



+gV A (r)TY {0 f r •ES + S t r-EO} + £-y^Tr {0 f r • E + f 0r • e} 



4 



flU 



where P is the momentum associated to the centre-of-mass coordinate. In Eq. ( |4"5[ ) the 
1/m corrections to Va, Vb and to pure gluonic operators as well as the higher order terms 
in the multipole expansion are not displayed. To the leading order the singlet and octet 
operators S, O are represented by the appropriate products of nonrelativistic heavy quark 
and antiquark spinors. The matching of S and O operators with the NRQCD spinors was 
done in |TIJ up to three loops for both the potentials and the normalization factors in 



OPE. In this lagrangian the singlet and octet potentials V s (r) and V (r) are treated as the 
corresponding Wilson coefficients in front of bilinear forms in S and O to the leading order in 
1/m. In ref. |10| the authors shown that this definition of static quark potential is consistent 



with the definition in terms of Wilson loop ([!]) . 

The other result of pNRQCD is the cancellation of renormalon ambiguity in the sum 
of heavy quark pole masses and the potential up to two loops, which is a confirmation of 
general consideration in QCD, that was first derived in |[42| . 

A new feature appears by the consideration of three-loop leading log matching of V and 
MS schemes. So, for the distance-dependent running coupling the result reads off 

( aMs(lA 2 ) 



a v (l/r ,fi) = Q!Ms(l/r ) { 1 + [a x + 2>y E (3 
+ 



V 



lE (4oiA, + 2/?i) + + 4 7 ^ 1 (3 2 + a 2 



47T 



16 7T 2 



(46) 



12 7T 



where the two-loop contribution was taken from [|7],[8| and it is coincides with (|I4]), of course. 
However, the three-loop term leads to the explicit dependence on the scale in the perturba- 
tive pNRQCD calculations, which has to be expected from the general note on the infrared 
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singularity observed by Appelquist, Dine and Muzinich ||12|| , that was rederived in pNRQCD 



by supplementing a certain infrared subtraction. This dependence was considered in [H| 
for two cases, when the scales of confinement Aqcd and binding energy mv 2 have the ar- 
rangements a) Aqcd 3> m v 2 or b) mv 2 ^> Aqcd- If «), the singlet potential of static quarks 
suffers from the nonperturbative effects, and it can be treated only after introduction of 
some model dependent terms coming from the ultrasoft gluons, which form the gluon sea 
in the heavy quarkonium, so that the sea has its excitations, and the characteristic excita- 
tion energy of gluelumps should replace the scale /x, that results in the scale-independent 
nonperturbative potential^. If b), the potential is purely perturbative. However, calculating 
the physical quantities such as the masses of bound states, we have to take into account 
the contributions coming from the perturbative ultrasoft gluons with the virtualities less 
than /i, which can produce a /i-dependent shift of energy, that should be cancelled with the 
/i-dependence in the potential fl46[) and, probably, in the heavy quark masses. In both cases, 
the perturbative calculations of singlet potential^ explicitly indicate the necessity to take 
into account the gluon degrees of freedom inside the heavy quarkonium. As was noted in 
]I(J , apparently, this feature is characteristic for the nonabelian theory (see the factor of Ca 
in front of log term in ( fl6|) ) . 

To our opinion, this dependence of potential on the ultrasoft gluon fields (the infrared 
singularity in terms of Appelquist, Dine and Muzinich) inside the heavy quarkonium nat- 
urally indicates the formation of gluon string between the heavy quarks at long distances. 
Indeed, expression (45) was derived under the following arrangement of scales: r ~ mv, 
mv 2 < fi < mv. So, if we put 

u 

fX = cr r, 

r 

with 

u 

v < u < 1, and a <C 



m z v 



then, perturbatively expanding in the small parameter ar, we get the linear correction to 
the potential in pNRQCD, so that 



where 



AVpnrqcd = Ak-r, 



so that we have dropped the scale dependence of strong coupling constant, since it is beyond 
the accuracy under study, and we have substituted the coulomb relation for the quark 
velocity inside the bound state u ~ CfQ-ms to the given order. Numerically, for the charmed 
quarks this perturbative contribution could be of the order of Ak ~ 0.1 GeV 2 . Thus, we can 



* Possible non-potential terms are discussed in [[To,]. 

TWe do not concern for the octet potential of static quarks in the present consideration, though 
some qualitative conclusions could be straightforwardly generalized from the siglet state to the 
octet one. 
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motivate the relation between the nonperturbative string and the three-loop scale dependent 
term in the pNRQCD potential. 

Indifferently of the arrangement for the confinement and binding energy scales, the in- 
troduction of such the string should remove the explicit dependence of full potential on 
the scale. This has been done above by introduction of unified (3 function of coupling in 
the V scheme. This solution of the problem qualitatively agrees with the consideration in 
pNRQCD, since, first, in the perturbative regime the contribution of log term is negligibly 
small as we see for the linear confining term of potential at short distances, and, second, 
at long distances the nonperturbative confining term is essential, where the string tension 
is the natural physical scale. In the static potential of QCD given above we do not con- 
sider possible 'nontrivial' excitations with the broken string geometry, where the break point 
moves on the string with the speed of light. Such the excitations would correspond to hybrid 
states with the gluelumps. Thus, we find that the QCD potential of static quarks in the 
form offered in the present paper has no conflicts with the current status of pNRQCD. 

However, to our opinion the problem can be more deep. The static potential, introduced 
by the Wilson loop, is renormalization group invariant, and it does not contain any separation 
between the potential gluons and the ultrasoft gluons forming the sea, since it gives the 
total energy of dynamical fields. In contrast, the pNRQCD introduced the singlet potential 
as the Wilson coefficient in front of four quark operator, so that it intrinsically operates 
with the separation of potential and sea, as well as the nonrelativistic quarks, which act as 
sources, so that some gluons with virtualities greater than /i are considered as emitted, while 
others with virtualities less than \x are included into the origin of sources, and the gluons 
with virtualities about mv mediate the potential interaction. Generally, this separation of 
heavy quarks, potential gluons and sea gluons in the operator product expansion can involve 
nonzero anomalous dimensions for the singlet pNRQCD-potential, say. This fact does not 
contradict with the OPE basis, but it reflects the point that the static potential of Wilson 
loop generally differs from the pNRQCD-potential. In addition, the ultrasoft gluon sea 
introduced in pNRQCD in terms of multipole interaction with local external chromoelectric 
and chromomagnetic fields is not a local object, indeed. 

A note should be done on the linear confining term of potential. In |TIJ a model of 
infrared behaviour was used, so that at long distances between the heavy quarks the ultrasoft 
correction was derived in the form of constant energy shift 5Vq and quadratic term a^T 2 . 
The corresponding conclusion was drawn to stress that the linear term could appear in a 
more complicated case of infrared behaviour. We show in the previous sections how this 
confinement regime can be reached. 

Recently, several papers were devoted to the calculations of ground states in the 

heavy quarkonia in the way, combining the pNRQCD potential with the nonperturbative 
corrections to the binding energy as they produced by the multipole expansion of QCD 
If| in the form of pNRQCD explicitly shown in (|4"5|). Ref. jPJ does not strictly estimate 



the gluon condensate effects in the multipole expansion, and it presents purely perturbative 
results. It follows the perturbative ground state mass technique as a mass definition that 
leads to the cancellation of the u = 1/2 renormalon that was considered in the approach 
of upsilon expansion introduced by Hoang et al. in [^6[ . So, in |B| the perturbative mass 
of B c meson was calculated on the base of perturbative expansion for the static potential 
with the leading approximation in the form of coulomb wave functions. As we see above 
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the perturbative potential suffers from the renormalon ambiguity. In order to remove this 
dependence on the choice of scale \i in the potential, the authors of [fH]] calculated the 
masses of J/^> and T in the same technique at the same point \i and inverted the problem 
on the heavy quark masses by equalizing the perturbative masses of ground states in the 
charmonium and bottomonium to the measured values. This procedure leads to the \x- 
dependent pole masses of heavy quarks as expressed by the series in a s {p)- We expect that 
such the procedure could cancel the renormalon with the accuracy about 50 MeV in the 
mass of hadron. As a results, the perturbative mass of B c has quite a stable value 



M pevt (B c 



6326±| 9 MeV, 



(47) 



in the range of 1.2 < fi < 2.0 GeV, which should be compared with the results in Table [D 1 
and the range of \x described above in the study of matching the perturbative potential with 
the full QCD one. The authors of [33| did not present the /i-dependent heavy quark masses. 
Nevertheless, due to the almost coinciding estimates of B c mass in (f47|) and Table |TJ, we 
expect that this dependence should be given by the form of 8V(fi). 

In refs. |45| the same technique for the perturbative contribution with the account for 
both the gluon condensate corrections in the multipole expansion of QCD and a small 
a B s log ot s term, was used to extract the heavy quark masses. The authors determined the 
'pole' mass, which is scale dependent, indeed, by putting fi = Cpa s mQ in the potential. As 
we understand, they introduced the mass suffered from the renormalon and got 

m b = 5022 ± 58 MeV, 

which is greater than we determine in the current presentation. The reason is quite evident. 
It is the energy shift SV(fi). The running MS mass quoted in j|5j is about 260 MeV greater 
than we find in the same order in a s for the relation between the pole and running masses. 



The difference becomes unessential by using the three-loop matching of the masses in [45 



however, the same correction will also decrease the value obtained in the spectroscopy with 



the full QCD potential. Thus, to our opinion the values of heavy quark masses given in f45 
should be kept with a large care. 



Finally, in [H7[ the dependence of potential on the finite heavy quark masses was con- 
sidered. This dependence is due to the smooth variation of number of active flavors in the 
expressions for the coefficients of perturbative (3 function as well as in the matching coeffi- 
cients of a v . As we have described above we use the step-like change of active flavor number, 
which infers implicit model-dependence, which is practically unavoidable in the case under 
study. 

As for the lattice simulations in QCD for the relevant problem, the review can be found 
in ref. We emphasize only that the lattice potential of static quarks is close to what 
is given by the Cornell model. A modern review of phenomenological potential models can 

1491 . The finite mass effects in the nonrelativistic bound states was 



be found in lectures 

recently considered at nest-to-leading order in |5(| and [51]. A next-to-next-to- leading order 



analysis of light quark mass effects in the heavy nonrelativistic QQ systems was given in 



Some applications of pNRQCD to the heavy quarkonia were done in [53 . 
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C. Leptonic constants 



In the static approximation for the heavy quarks the calculation of leptonic constants for 
the heavy quarkonia with the two-loop accuracy involves the matching of leptonic currents 
in NRQCD with the currents of full QCD, 

JT D = Ql,Q, 3 NRQCD = xV; l 0) 

with the relativistic quark fields Q and their nonrelativistic two-component limits of anti- 
quark x an d quark <f>, = a u — v v {a ■ v), and v is the four- velocity of heavy quarkonium, 
so that 

jQCD _ )C(fJ, hgld ; //f act ) ■ jNRQCD ^ 

where the scale //hard determines the normalization point for the matching of NRQCD with 
full QCD, while //f act refers to the point of perturbative calculations in NRQCD. Using 
the matching of potential for the static quarks in QCD with the two-loop perturbative 
potential, we argue that the most appropriate choice of scale relevant to the charmonium 
and bottomonium is 



//fact = //soft = 1.3-2 GeV. (48) 
For the heavy quarkonium composed by quarks of the same flavour the Wilson coefficient 



/C is known up to the two-loop accuracy [B4LR3lp5l56 



w, \ -I 8a^ S (/i h ard) /«^ S (/ihard)\ , ,. 

A- (//hard; //fact) = 1 ~ ~ V C 2 (//hard 5 //fact J , (49) 

on \ 7T / 



and C2 is explicitly given in The additional problem is the convergency of ( [491 ) at the 

fixed choice of scales. So, putting //hard = (1 — 2)m& and (fHf ) we find a good convergency of 
QCD corrections for the bottomonium and estimate its leptonic constant defined by 

(0\J? CD \r,\) = e x JyM r , 

where A denotes the polarization of vector state e u , so that 

fr = 685 ± 30 MeV, 

while the experimental value is equal to /£ XP = 690 ± 13 MeV 0. 

As we can see in Fig. [7] the variation of hard scale in broad limits leads to existence of 
stable point, where the result is slowly sensitive to such the variation. The stability occurs 
at //soft ~ 2.6 GeV, where the perturbative potential is still close to the potential of static 
quarks at the distances characteristic for the lS'-level of bb. 

The estimate of leptonic constant for the charmonium J/ip is more sensitive to the 
choice of factorization scale. Indeed, the size of this system, (rgcfis)) ~ 0.42 fm, makes more 
strict constraints on // fact « 1.3 — 1.5 GeV, since at higher scales the perturbative potential 
significantly deviates from the potential of static quarks in QCD in the region of bound cc 
states, while at lower scales the perturbative potential in two loops does not match the QCD 
potential in all of the form. 
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fr, GeV 




FIG. 7. The value of leptonic constant for the vector ground state of bottomonium versus the soft 
scale. The dashed line represents the choice of //hard = 2m&, while the solid line does //hard = m b- 
The horizontal shaded band gives the experimental limits. 
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FIG. 8. The value of leptonic constant for the vector ground state of charmonium versus the 
soft scale. Shaded region is restricted by the dashed line representing the choice of //hard — 

1.07m,, 

and the solid line with //hard = 0.93m c . The horizontal shaded band gives the experimental limits. 
The additional curves give //hard = 1.26m c (dashed line) and //hard = 0.87m c (solid line). 
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Another problem is the energy shift SV(fi) = 1.0—1.2 GeV, which essentially renormalizes 
the pole mass of charmed quark: m^ olc = 1.968 — 2.068 GeV. This shift does not perturb 
the mass of the ground state, but it is significant for the value of wave function at the 



origin. So, following the well-adjusted scaling relation for the leptonic constants |57|| , we put 
P(/i) = k\I/(0) mP° le (/i)/m c and use it in the calculations of the leptonic constantfj]. We get 

f m = 400 ± 35 MeV, 

to compare with the experimental value fjj^ = 409 ± 15 MeV. 

In Fig. [| we see that again the stability point can be reached in the variation of /ihadr at 
reasonable value of /i so f t ~ 1.35 GeV. However, the stability takes place in the narrow region 
of /ihadr close to the charm quark mass. 

At present, the matching condition for the heavy quarkonium composed by the quarks 
of different flavors, be, is known to one loop, only [fT7| , |5T| . So, for the pseudoscalar state we 
have 



^(/Uiard! /^fact) 



af S (//hard) 
7T 



m b -m c m b 
in — 



mi, + m c 



(50) 



which is independent of the factorization scale. The matching of perturbative potential to 
the one-loop accuracy with the QCD potential of static quarks at r ~ 0.3 — 0.4 fm relevant 
to the ground state of B c meson [j59|1 , is rather questionable, since the deviation in the forms 
of potentials is quite sizable. In addition we have to pose /tf act = /ihard, because we cannot 
distinguish these scales, while the nonzero anomalous dimension to two loops is not taken 
into account. Nevertheless, we can put /ihard = 1.3 — 1.8 GeV and neglect 5V, which is 
beyond the actual control in the one-loop accuracy. Indeed, as we see in Fig. |5] the one loop 
value of energy shift for the matching of perturbative and QCD potentials is quite small 
at the large virtualities about 2 GeV, and it can be neglected, while at smaller virtualities 
the form of perturbative potential is close to that of given by QCD only in the short range 
of distances r = 0.1 — 0.25 fm, hence, the results on the matching are not reliable for the 
extracting the heavy quark masses from the parameters of bound states. So, we estimate 

f Bc = 400 ± 45 MeV, 

to compare with the estimates in the SR, where /f* = 400 ± 25 MeV |7|,|0j. 



TABLE IV. The ratios of leptonic constants for the heavy quarkonia as predicted in the present 
paper (K 2 0) in comparison with the experimental data. 





QCD, K 2 


exp. 


fr(ns)/ fr 


QCD, K 2 


exp. 


2S 


0.55 


0.48 ± 0.07 


2S 


0.47 


0.47 ± 0.03 


3S 


0.32 


0.25 ± 0.06 


3S 


0.34 


0.36 ± 0.02 



£ Solving the Schrodinger equation with the shifted masses and potential, we check that this mass 
dependence of wave function is valid with the accuracy better than 6%, so we put k = 0.95. 
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Finally, we present the ratios of leptonic constants for the excited nS'-levels of bb and cc 
in Table |I^ in comparison with the experimental data. We see that the predictions are in a 
good agreement with the measured values. For completeness, we predict also the constant 
of 2S'-level in the be system 

f Bc (2S) = 280 ± 50 MeV, 
which agrees with the scaling relation p7 . 



Thus, we have analyzed the estimates following from the potential of static quarks in 
QCD for the masses of quarks and heavy quarkonia as well as for the leptonic constants, and 
found both the good agreement with the experimental data available and the consistency 
with the QCD sum rules. 

IV. CONCLUSION 

We have derived the potential of static heavy quarks in QCD on the base of known limits 
at short and long distances: the asymptotic freedom to the three loop accuracy and the 
confinement regime. The inputs of potential are the coefficients of perturbative (3 function, 
the matching of MS scheme with the V scheme of potential, the normalization of running 
coupling constant of QCD at /i 2 = m| and the slope of Regge trajectories, determining the 
linear term in the potential. Thus, the approach by Buchmiiller and Tye has been modified 
in accordance with the current status of perturbative calculations. 

In the static limit the two-loop improvement of coulomb potential results in the signif- 
icant correction to the j3 function for the effective charge, A/3/ (3 ~ 10% as shown in Fig. 
|3|. This correction is important for the determination of critical values of charge, i.e. the 
value in the intermediate region between the perturbative and nonperturbative regimes. 
Moreover, the two-loop matching condition and the three-loop running of coupling constant 
normalized by the data at the high energy of uiz determine the region of energetic scale for 
changing the regimes mentioned above. This scale strongly correlates with the data on the 
mass spectra of heavy quarkonia. So, it is connected with the splitting of masses between 
the IS* and 25' levels. We stress that the consistent consideration of two- loop improvement 
gives the appropriate value of effective coulomb coupling constant as it was fitted in the 
Cornell model of potential. This is achieved in the present paper in contrast to the one- 
loop consideration by Buchmiiller and Tye, who found the value of Aqcd inconsistent with 
the current normalization at high energies. So, the two-loop improvement gives the correct 
normalization of effective coulomb exchange at the distances characteristic for the average 
separation between the heavy quarks inside the heavy quarkonium and determines the de- 
viations at short distances r < 0.08 fm (see Fig. [|), that is important in the calculations of 
leptonic constants related with the wave functions at the origin. 

Other corrections to the potential of heavy quarks are connected with the finite mass 
effects and cannot be treated in the framework of static approximation. For example, the 
spin-dependent forces, relativistic corrections and specific non-abelian potential terms^] in 
the heavy quarkonium should be taken in the analysis of mass spectra. A magnitude of 



5 They have the form of ot 2 s /r 2 with the factor given by the inverse heavy quark masses. 
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leading nonstatic corrections can be evaluated by the characteristic shifts of levels due to 
the hyper-fine splitting of >S-wave levels in the heavy quarkoniaQ- So, we conservatively 
evaluate the uncertainty of heavy quark mass analysis 8m ~ 80 MeV. 

Thus, the non-abelian term of potential ce 2 /r 2 , say, has the factors in the form of I/tuq, 
and it is equal to zero in the static limit itlq — > oo, while the uncertainty in the heavy 
quark masses due to the omission of such the terms is estimated in the paragraph above. 
Formally, if we consider the perturbation theory for the calculation of bound state levels in 
the heavy quarkonium with the coulomb functions taken as the leading approximation, which 
is not a scope of our consideration, then the mentioned non-abelian potential contributes 
in the same order in a s as the two-loop corrections to the matching of perturbative static 
potential ~ af, since the averaging of 1/r 2 results in a 2 mQ factor. However, the two-loop 
effects are important for the consistent consideration of static potential and the high energy 
normalization, i.e. these corrections are significant in the running of effective charge in the 
potential from the high energies to the scale relevant to the heavy quark bound states even 
in the static limit, while the nonstatic contributions can be consistently neglected in the 
numerical analysis. We see that our consideration is consistent in the static approximation, 
which we have addressed in the present paper. 

The matching of two loop perturbative potential with the QCD potential of static quarks 
has been performed to get estimates of heavy quark masses, which can be compared with 
the results of QCD sum rules. A good agreement between two approaches has been found. 

The recent determinations of heavy quark masses in Refs. [[3^-|34| were done in the 
framework of QCD sum rules, which is a systematic approach, indeed. It is based on 
the separation of short- distance region from the nonperturbative effects at some values of 
parameters defining the scheme of calculations in the sum rules. In this approach the 
nonperturbative terms are given in the form of quark-gluon condensates contributing with 
corresponding short- distance Wilson coefficients, so that as was shown in |37|, a numerical 
contribution of gluon condensate term in the sum rules is negligibly small in comparison 
with the perturbative part. However, it would be incorrect to think that these explicit 
contributions suppressed in some region of parameters are the only terms caused by the 
nonperturbative infrared dynamics of QCD. Indeed, neglecting the condensate terms, we 
find that the perturbative correlators suffer from the renormalon ambiguity, which implies 
that the perturbative expansion in series of a s is asymptotic, and the summation of series 
depends on a method used. The physical reason for such the divergency and ambiguity 
is the infrared singularity in the QCD coupling constant. This singularity is regularized 
by introducing the threshold mass parameters free of renormalon. Such the approach is 
independent of any assumptions on the gluon condensate, since, generally the pole mass 
renormalon and the gluon condensates are different issues. 

So, the perturbative pole mass used in the QCD sum rules is not well defined quantity, 
and some relevant quantities are introduced in Refs. |32|-|3~4}| . These quantities are con- 
structed from the perturbative pole mass of heavy quark with specific infrared subtractions, 
which are treated independently of the quark-gluon condensates. These constructions are 
author-dependent. Though the authors of subtracted masses gave some physical motiva- 



**The splitting is about 100 MeV or less. 
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tions, which are more or less strict, but justified. These infrared subtractions imply the 
introduction of infrared regulators. 

In the present paper the unified (3 function for the effective charge in the potential is 
considered, and its definition supposes the infrared stability. Thus, we see that the analysis 
of heavy quark masses in both the QCD sum rules and potential approach involves the 
consideration of relevant effects caused by the infrared dynamics of QCD, though the explicit 
constructive procedures are certainly different, but they have similar inherent uncertainties. 

The calculated mass spectra of heavy quarkonia and the leptonic constants of vector 
nS-levels are in agreement with the measured values. The characteristics of B c meson have 
been predicted. 
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